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Abstract. CMB data analysis is in general done through two main steps : map-making of the time 
data streams and power spectrum extraction from the maps. The latter basically consists in 
the separation between the variance of the CMB and that of the noise in the map. Noise 
must therefore be deeply understood so that the estimation of CMB variance (the power 
spectrum) is unbiased. I present in this article general techniques to make maps from time 
streams and to extract the power spectrum from them. We will see that exact, maximum 
likelihood solutions are in general too slow and hard to deal with to be used in modern 
experiments such as Archeops and should be replaced by approximate, iterative or Monte- 
Carlo approaches that lead to similar precision. © 2003 Academie des sciences 

Cosmology / Cosmic Microwave Background / Data Analysis 



1. Introduction 

The cosmological information contained in the Cosmic Microwave Bacground (CMB) anisotropies is 
encoded in the angular size distribution of the anisotropies, hence in the angular power spectrum and noted 
Ci. It is of great importance to be able to compute the Ci spectrum in an unbiased way. The simplest 
procedure to obtain the power spectrum is to first construct a map of the CMBA from the data timelines 
giving the measured temperature in one direction of the sky following a given scanning strategy on the sky, 
this is known as the map-making process ; then extract the Ci from this map, this is the power spectrum 
extraction. Various effects usually present in the CMB data make these two operations non trivial. The 
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major effect being related to the unavoidable presence of instrumental and photon noise. Noise in the 
timelines is correlated and appears as low frequency drifts that are still present in the map. A good map- 
making process minimizes these drifts, but in most cases, they are still present in the map. They have to 
be accounted for in the power spectrum estimation as the signal power spectrum is nothing but an excess 
variance in the map at certain angular scales compared to the variance expected from the noise. The CMBA 
power spectrum will therefore be unbiased only if the noise properties are known precisely. 

This article presents the usual techniques that allow an unbiased determination of both the CMBA maps 
and power spectrum. In Sect. |3 we will describe the data model and the data statistical properties required 
for the techniques presented here to be valid. Sections l3~l and l4l respectively deal with map-making and 
power spectrum estimation techniques. 

2. Data model 

The initial data are time ordered information (TOI) taken along the scanning strategy pattern of the 
experiment. The detector measures the temperature of the sky in a given direction through an instrumental 
beam. This is equivalent to say that the underlying sky is convolved with this instrumental beam and that 
the instrument measures the temperature in a single direction of a N p pixellised convolved sky noted T. 
The Nt elements TOI noted d may therefore be modelled as: 

d = ,4 T + n (1) 

The pointing matrix A relates each time sample to the corresponding pixel in the sky. A is a N t x N p matrix 
that contains a single 1 in each line as each time sample is sensitive to only one pixel is the convolved sky 1 . 
The noise TOI n in general has a non diagonal covariance matrix N given by 2 : 

N=(n- n*) (2) 

The most important property of the noise, that will be used widely later is that it has to be Gaussian and 
piece-wise stationary. Both assumptions are crucial as they allow major simplifications of the map-making 
and power spectrum estimation problems, namely Gaussianity means that all the statistical information on 
the noise is contained in its covariance matrix and stationarity means that all information is also contained 
in its Fourier power spectrum leading to major simplifications of the covariance matrix : the noise depends 
only on the time difference between two samples and N is therefore a Toeplitz matrix iVy = Nu_ji 
completely defined by its first line and is very close to be circulant 3 . Such a matrix is diagonal in Fourier 
space. Its first line is given by the autocorrelation function of the noise, that is the inverse Fourier transform 
of its Fourier power spectrum (* is the convolution operator) 4 : 

(3) 



N i0 = {n*n) = (F- 1 \T(n 



3. Map-making techniques 

The map-making problem is that of finding the best estimate T of T from Eq.^ given d and A. The 
noise n is of course unknown. We will address the two main approaches to this problem, the first being 
the simplest one and the second one being the optimal one. An excellent detailed review on map-making 
techniques for the experts is | Stompor, 2002 1. 



'Different forms for A can however be used in case of differential measurements or more complex scanning strategies. 
2 the symbols { ) mean that we take the ensemble average over an infinite number of realisations. 
3 Saying that the matrix is circulant is an additionnal hypothesis, but a very good approximation for large matrices. 
4 JF denotes the Fourier transform (in practice, a FFT algorithm is used). 
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3.1. Simplest map-making : coaddition 

The simplest map-making that one can think about is to neglect the effects of the correlation of the noise. 
One can just average the data falling into each pixel without weighting them. This procedure is optimal (it 
maximises the likelihood) if the noise in each data sample is independant, that is, if the noise is white. In a 
matrix notation, this simple map-making can be written: 

T = [A 1 ■ A}' 1 ■ A 1 ■ d (4) 

where the operator A t just projects the data into the correct pixel and \A l ■ A] counts the sample falling into 
each pixel. This simple map-making has the great advantage of the simplicity. It is fast (oc N t ) and robust. 

However, in the case of realistic correlated noise, the low frequency drifts in the timelines induce stripes 
in the maps along the scans of the experiment. These stripes are often much larger than the CMBA signal 
that is searched for and therefore should be avoided. Various destriping techniques have been proposed to 
avoid these stripes. A method exploiting the redundancies of the Planck mission 5 scanning strategy has 
been proposed by | Delabrouille, 1998 1 and extended to polarisation by |R evenu et al., 2000| . This kind 
of method aims at suppressing the low frequency signal by requiring that all measurements done in the 
same direction at different instant coincide to a same temperature signal. Another method has recently been 
proposed for the Archeops 6 data analysis and estimates the low frequency drifts by minimizing the cross- 
scan variations in the map due to the drifts | Bourrachot et al., 2003 1. The simplest method for removing 
the low frequency drifts before applying simple map-making is certainly to filter the timelines so that the 
resulting timeline has almost white noise. The filtering can consist in prewhitening the noise or directly 
setting to zero contaminated frequencies. The computing time (CPU) scaling of the filtering + coaddition 
process is modest and dominated by filtering (oc Nt log Nt). This method however removes also part of the 
signal on the sky and induces ringing around bright sources which has to be accounted for in later processes. 

3.2. Optimal map-making 

The most general solution to the map-making problem is obtained by maximizing the likelihood of the 
data given a noise model |Wright, 1996||Tegmark, 19971 . As the noise is Gaussian, its probability distribu- 
tion is given by the Nt dimensionnal Gaussian: 



P(n) = -pr exp 

| (2tt)^JV 1 1/2 



1 i i 
-n* • N n 

2 



(5) 



Assuming no prior on the sky temperature, one gets from Eq.^the probability of the sky given the data: 

(6) 



P(Tld) oc P(dlT) oc --r exp 

|(27r)^A| 1/2 



-~(d - A • T)* • A" 1 • (d - A • T) 



Maximizing this probability with respect to the map leads to solving the linear equation: 

A* • A" 1 • A ■ T = A* • A" 1 • d (7) 

with solution 7 : 

t = {A t ■ N- 1 -Ay 1 ■ A* ■ N' 1 ■ d (8) 

One therefore just has to apply this linear operator to the data timeline to get the best estimator of T, note 
that T is also the minimum variance estimate of the map. The covariance matrix of the map is: 

N = (A* ■ N- 1 ■ A)' 1 (9) 



http://astro.estec.esa. nl/Planck71 



http : / /www . archeops .org/ 



One can remark here that simple map-making is equivalent to optimal map-making if the noise covariance matrix is diagonal, 
which is consistent to what was said before. 



3 



J.-Ch. Hamilton 



THE COSMIC MICROWAVE BACKGROUND 



Problems arise when trying to implement this simple procedure, the timeline data and the maps are in 
general very large : the typical dimensions of the problem are N t ~ 6 x 10 7 and N p ~ 10 5 for Archeops. 

The maximum likelihood solution requires both A -1 and (A* • TV" 1 • A) 1 which are not easy to de- 
termine. Two approaches can be used at this point: one can try to make a brute force inversion of the 
problem, relying on huge parallel computers or one can try to iteratively approach the solution, hoping that 
convergence can be reached within reasonnable time. 

3.3. Brute force inversion 

The brute force optimal map-making parallel implementation is freely available as the MADCAP [ Borrill, 1 999 1 
package. It is a general software designed to produce an optimal map for any experiment by solving directly 
Eq.Q The use of this package requires the access to large parallel computers. 

The only assumption that is done in MADCAP map-making is that the inverse time-time noise covariance 
matrix can be obtained directly without inversion from the noise Fourier power spectrum: 



\H*)\ 2 



(10) 



This assumption is not perfectly correct on the edges of the matrix but leads to a good estimate of the inverse 
time covariance matrix for the sizes we deal with. This allows this step to scale as N t \ogN t operations 
rather than the Nf required by a Toeplitz matrix inversion. In most cases, the time correlation N T length is 
less than the whole timestream N t so that N is band-diagonal. For Archeops, we have N T ~ 10 4 . 

The next step is to compute the inverse pixel noise covariance matrix A/" -1 = (A* ■ N^ 1 ■ A) and the 
noise weighted map A t ■ N^ 1 ■ d, both operations scale as Nt x N T when exploiting the structure of A and 
N. The last step is to invert JV^ 1 and multiply it by A 1 ■ N^ 1 ■ d to get the optimal map. Unfortunately, 
A/" -1 has no particular structure that can be exploited and this last step scales as a usual matrix inversion 
oc N^ and largely dominates the CPU required by MADCAP for the usual large datasets (eg. Archeops). 

We can remark here that MADCAP provides the map covariance matrix J\f for free as a byproduct. This 
matrix is crucial for estimating the power spectrum as will be seen in sectionl4~l 

3.4. Iterative solutions 

The other possibility is to solve Eq.Q through an iterative process such as the Jacobi iterator, or more 
efficiently a conjugate-gradient | Press et al., 1988 1. Both converge to the maximum likelihood solution. 

The use of the Jacobi iterator for solving for the maximum likelihood map in CMB analysis was first 
proposed by | Prunet et al., 20 00 1. The basic algorithm is the following. We have to solve the following 
linear system (see Eq.0: 

r-x = y (11) 
The Jacobi iterator starts with an approximation Ao of and iterates to improve the residuals R: 

A -T = I-R (12) 

In order to converge, the algorithm requires the first approximation to be good enough so that the eigenval- 
ues of R are all smaller than 1 (a good estimate in general is Ao = [diagT] -1 ). We can therefore expand: 

r- 1 = (/ - Ry 1 • A = (7 + R + R 2 + •••)• A (13) 

Lets us define A„ = (/ + R + R 2 + ■ ■ ■ + R n ) ■ A so that = lim„_ >00 A„. We have the relationship 
Aj_|_i = Aq + R ■ Aj. If we define xj = Aj • y, it is straightforward to show that: 



- XJ = A • (y - r • Xj ) (14) 
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which defines the Jacobi iterator. When going back to the usual CMB notation for maps and timelines, one 
gets: 

T i+ i - Tj = [diag (A* • A/" 1 • A)] ~ 1 ■ A* ■ N^ 1 • (d - A ■ Tj) (15) 

which looks rather complicated but is in fact very simple to implement: the operation A ■ Tj just consists 
in reading the map at iteration j with the scanning strategy (oc N t ), and the matrix diag (A* ■ N^ 1 ■ A) is 
just the white noise level variance divided by the number of hits in each pixel. It is diagonal and therefore 
does not require proper inversion. The only tricky part here is the multiplication TV -1 • (d — A ■ Tj) given 
the fact that N^ 1 is unknown. As the noise is stationary, N is Toeplitz and circulant 8 , the multiplication by 
TV -1 can be done in Fourier space directly through: 



iV" 1 x~ T 



-l 



l^(x)| 2 



(16) 



which requires Nt log Nt operations. Finally, each iteration is largely dominated by the latter so that the 
final CPU time scales like Nu x Nt log Nt where Nu is the number of iterations. 

Unfortunately the convergence of such an iterator is very slow and makes it rather unefficient as it is. A 
significant improvement was proposed by |Dore et al., 2001 1 in the publicly available software MAPCUMBA. 
They noted that the convergence was actually very fast on small scales (compared to the pixel) but that the 
larger scales were converging slowly. They proposed a multigrid method where the pixel size changes at 
each iteration so that the global convergence is greatly accelerated (see Fig. 7 of |Dore et al., 2001 1), making 
this iterative map-making really efficient. A conjugate gradient solver instead of the Jacobi iterator is imple- 
mented in the software Mirage |Yvon et al., in prep!) and accelerates again the convergence significantly. A 
new version of MAPCUMBA also uses a conjugate gradient solver, as well as MADmap |Cantalupo, 2002 1. 

If obtaining an optimal map is now quite an easy task using an iterative implementation (the presence 
of strong sources, such as the galactic signal however complicates this simple picture), they do not provide 
the map noise covariance matrix M — (A* • N^ 1 • A) which is of great importance when computing 
the CMB power spectrum in the map in order to be able to make the difference between noise fluctuations 
and real signal fluctuations. The only way to obtain this covariance matrix using these iterative methods is 
through large Monte-Carlo simulation that would reduce the advantage of iterative map-making compared 
to brute-force map-making. 

3.5. map-making comparisons 

The precision of the MADCAP, MAPCUMBA and Mirage implementations are shown in Fig with 
the same CMB and noise simulation based on Archeops realistic conditions. The three resulting maps 
were kindly provided by | Filliatre, 2002 1. The six maps on the left are respectively from top left to bottom 
right : initial CMB fluctuation, coaddition of the timeline without filtering, coaddition of the timeline with 
white noise only (ie the true optimal map that has to be reconstructed), MADCAP residual map (difference 
between MADCAP reconstructed map and the white noise map), MAPCUMBA residual map and Mirage 
residual map. All maps are shown with the same color scale. The first remark that can be done is that the 
stripes are indeed a real problem and that straight coaddition is not to be performed. The three different 
optimal map-making codes give very similar results, especially MADCAP and MAPCUMBA. In all cases, 
as can be also seen in the right panel of Fig. [2 the residuals are much smaller than the CMB fluctuations 
that are searched for. The three map-making implementation can therefore be considered are unbiased 9 . 

Finally one can summarize the comparison as following, iterative and brute-force optimal map-making 
give very similar results as far the optimal map is concerned. The brute force inversion provides the map 
noise covariance matrix for free which is a major point as will be seen in next section. The computer 



8 again, it is not exactely circulant but it is an excellent approximation as the matrix is large 

'Let's note that the noise model that was used for MADCAP is the true one, not an estimation. This makes however little difference. 
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requirements are however much larger than for iterative map-making. The latter should therefore be used 
when the power spectrum estimation can be carried out without the knowledge of the map noise covariance 
matrix, in general using a Monte-Carlo technique (see next section). In this case, one should seriously 
consider the filtering + coaddition map-making that is by far the fastest but removes part of the signal. This 
is however accounted for (see section l4~2l also using a Monte-Carlo technique. 




Fi gure 1 : The six maps on the right show a comparison of results from different map-making implementations on 
the same simulation (typical of Archeops data). All maps are in Healpix pixellisation [Gorski et al. 1998 ] and have 
the same color scale ranging from -500 /iK to 500 /j,K from dark blue to dark red, green corresponds to zero. The 
histograms of the residuals is shown on the right and is more than three times smaller than the actual CMB fluctuation. 



4. Power spectrum estimation techniques 

We know want to compute the power spectrum Cg of the map T whose noise covariance matrix J\f might 
be known or not depending on the method that was used before to produce the map. The map is composed 
of noise and signal (from now on n is the noise on the map pixels): 

T = s + n (17) 

The signal in pixel p can be expanded on the Yt m (0 P , <fi p ) spherical harmonics basis : 

oo £ 

8 p=/ j ^ at m BtYi m (6 p ,<f> p ) (18) 

1=0 m=-i 

where Bg stands for the beam 10 . If the CMBA are Gaussian, the variance of the ag m , called the angular 
power spectrum and denoted Cg contains all the cosmological information : 

(fllmflf'm') = Cg5u/S mm ' (19) 

The map covariance matrix (assuming no correlation between signal and noise) is: 

M = (T ■ T*) = (s • s*) + (n • n*) (20) 
= S + N (21) 

and the signal part is related to the Cf. 

S pp > = (s p s p <) = ^2 —£— c t B i p e(Xpp') (22) 



It is the Legendre transform of the instrumental beam under the assumption that it is symmetric. 
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where Xvv' = cos ( u p ' u p' )> u p being the unit vector towards pixel p and Pi are the Legendre polynomials. 

One therefore has a direct relation between the map and noise covariance matrices and the angular power 
spectrum: 

M = N + J2 '^-C i B 2 e P e ( Xpp > ) (23) 
t 

The power spectrum estimation consists in estimating Ce from T and N (that can be unknown) using this 
relation. 

4.1. Maximum likelihood solution 

Full details concerning this can be found in |Bond, Jaffe and Knox, 1998||Tegmark, 1997| . As for the 
map-making problem, the maximum likelihood solution proceeds by writing the probability for the map 
given its covariance matrix assuming Gaussian statistics 11 : 



P(Q|T)«P(T|C/) = (27r) 



-N p /2 



exp 



[(T* • AiT 1 -T) +Tr(lnM)] 



and we therefore want to maximize the likelihood function through = 0: 

L(C t ) =~ [(T* • Ar 1 • T) + Tr (InM)] 



Tedious calculations lead to the solution: 

Ct=J2 F U' xTr 

v 

where F is the Fisher matrix: 

F U j = Tr 



(T-T* -N) ■ M~ 



as 
ac. 



■ M~ 



as 

dC~i 



■ M~ 



OS 
dCh> 



■ M~ 



(24) 



(25) 



(26) 



(27) 



Eq. \26\ let Ci appear in both sides (in M) in an uncomfortable way and therefore cannot be solved sim- 
ply. The method usually used |Bond, Jaffe and Knox, 1998||Borrill, 19991 is the Newton-Raphson iterative 
scheme: One starts from an initial guess for the binned power spectrum 12 C^ ) and iterates until conver- 
gence following: 

C (i+1) = C w + <5C (28) 



with: 



SC = - 



' d 2 L 




-l 


dL 


dC 2 






dC 



(29) 



C=Ci 



the likelihood L being that of Eq.[25] Convergence is usually reached after a few iterations. The explicit 
form of the derivatives of Eq.l29lis: 



dL 
dC b 



dS 



8 2 L 



dc h dc h , 



dC b 
dS 



■ M~ x ■ m Tr 



M~ 



BCi 



■ M~ 



dS 
dC h , 



■ M -1 rn 



as 
ac b 



as as 
ac h ac h , 



(30) 
(31) 



"The trace appears from |M[ _1 = exp [— Tr (InM)] as the trace is invariant. 

12 Binned power spectrum means that we do not consider one single mode i but a bin in £ as we do not have access in general to all 
modes due to incomplete sky coverage. 
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where the index b denotes the bin number. 

Each iteration will then require a large number of large matrix operations forcing such an algorithm 
to be implemented on large memory parallel supercomputers. MADCAP [Borrill, 1999 1 is the common 
implementation of this algorithm and scales as 2(Nt, + §)iVp operations per iteration. The CPU/RAM/Disk 
problem is therefore even cruder for the power spectrum than for the map-making. This algorithm leads 
to the optimal solution accounting correctly for the noise covariance matrix and additionnaly provides the 
likelihood shape for each bin through the various iterations allowing to a direct estimate of the error bars. 

4.2. Frequentist approaches 

An alternative approach to power spectrum estimation is to compute the so called pseudo power spec- 
trum (harmonic transform of the map, noted C£) and to correct it so that it becomes a real power spec- 
trum. This approaches have been proposed and developped in | |Hivon et al., 2001] |Szapudi et al. 2 001 1. 
The harmonic transform of the map differs from the true Cg in various ways (we follow the notations 
from |Hivon et al., 2001 1): The observed sky is convolved by the beam and by the transfer function of the 
experiment so that the observed power spectrum is BfFgCg, where Bg characterizes the beam shape in har- 
monic space and Fg the filtering done to the data by the analysis process (that may also include electronic 
filtering by the instrument itself). The observed sky is in general incomplete (at least because of a Galactic 
cut) leading to the fact that the Cg measured are not independant as they are convolved in harmonic space by 
the window-function I White and Srednicki 19951 . We therefore have access to J2t< MggrB^,FgrCii where 
Mggi is the mode mixing matrix. Finally, the noise in the timelines projects on the sky and adds its contri- 
bution Ni to the sky angular power spectrum. At the end, the map angular power spectrum, the pseudo-Cg 
is related to the true Cg via: 

C e = J2 M u > B% Fgi C v + Ng (32) 

The frequentist methods propose to invert Eq.[32]making an extensive use of Monte-Carlo simulations 
(details can be found in |Hivon et al., 200 ll ): 

• the pseudo power spectrum of the map Cg is computed by transforming the map into spherical har- 
monics (generally using Healpix pixellisation and the anaf ast procedure available in the Healpix 
package IGorski et al. 19981 V 

• The mode mixing matrix is computed analytically through: 

^-^E^ + DwJ'o 1 "l) 03) 

where Wg is the power spectrum of the window of the experiment (in the simplest case 1 for the 
observed pixels and elsewhere, but more complex weighting schemes may be used, as in Archeops 
IBenoit et al., 2003| or WMAP | Hinsha wltal., 2003| ). In the SpICE approach [Szapudi et al. 200 1| , 
the Mggi inversion in harmonic space is replaced by a division in angular space which is mathemati- 
cally equivalent. 

• The beam transfer function is computed from a Gaussian approximation or the legendre transform of 
the beam maps or a more complex modelling if the beams are asymetric, such as in | Tristr arrTet al., in prep| 

• The filtering transfer function is computed using a signal only Monte-Carlo simulation (it should 
include the pre-processing applied to the time streams). Fake CMB sky are passed through the in- 
strumental and analysis process producing maps and pseudo power spectra. The transfer function is 
basically computed as the ratio of the input model to the recovered ensemble average. An important 
point at this step is to check that the transfer function is independent of the model assumed for the 
simulation. Let's also remark that using a transfer function that depends only on I is a bit daring 
as the filtering is done in the scan direction which, in general corresponds to a particular direction 
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on the sky. This approximation however seems to work well and has been successfully applied to 
Boomerang | Netterfi eld et al., 2002 1 and Archeops |B enoit et al., 2003*1 ). 

• The noise power spectrum is computed from noise only simulations passing again through the instru- 
mental and analysis process to produce noise only maps and pseudo power spectra. The noise power 
spectrum is estimated from the ensemble average of the various realisations. 

• Error bars are computed in a frequentist way by producing signal+noise simulations and analysing 
them as the real data. This allows to reconstruct the full likelihood shape for each power spectrum 
bin and the bin-bin covariance matrix. 

Such an approach based on simulations has the advantage of being fast: each realisation basically scales as 

3/2 

oc Nt log Nt for the noise simulation and map-making (if filtering + coaddition is used) and oc N p for the 
CMB sky simulation and pseudo power spectrum computation. An important advantage of such a method 
is the possibility to include in the simulation systematic effects (beam, pointing, atmosphere, ...) that would 
not be easily accountable for in a maximum likelihood approach. 

4.3. Cross-power spectra 

When several photometric channels are available from the experimental setup, it is possible to compute 
cross-power spectra between the channel rather than power spectra of individual channel or of the average 
of all channels. This has the advantage of suppressing the noise power spectrum (but not its variance of 
course) that is not correlated between channels and leaving the sky signal unchanged. The cross-power 
spectrum of channels i and j is defined as: 



The cross-power spectrum method can easily be associated with the frequentist approach simplifying sig- 
nificantly it implementation as one the most difficult part, the noise estimation, is now less crucial as noise 
disapears and cannot bias the power spectrum estimation. This has been successfuly applied in the WMAP 
analysis |Hins haw et al., 20031 . 

4.4. Which power spectrum estimator should be used ? 

The maximum likelihood approach is undoubtly the best method to use if possible, but its CPU/RAM/Disk 
requirements are such that in practice, with modern experiments, it is very difficult to implement. It should 
however be considered to check the results on data subsets small enough to make it possible. The frequentist 
approaches are much faster and provide comparable precision in terms of error bars and permit to account 
for systematic effects in a simple manner. The tricky part is however to estimate the noise statistical prop- 
erties precisely enough. The same difficulty exists however in the maximum likelihood approach where the 
noise covariance matrix has to be known precisely. It is generally directly computed in the map-making 
process from the time correlation function, thus displacing the difficulty elsewhere. In any case, the noise 
model has to be unbiased as the final power spectrum is essentially the subtraction between the pseudo 
power spectrum of the map and the noise power spectrum. Estimating the noise properties is a complex 
problem mainly due to signal contamination and pixellisation effects. A general method for estimating 
the noise in CMB experiment is proposed in | Amblard & Hamilton, 2003 1 and was successfuly applied for 
the Archeops analysis [Beno rt et al., 200 3). When multiple channels are available, the frequentist approach 
applied on cross-power spectra is certainly the simplest and most powerful power spectrum estimation 
technique available today as it reduces the importance of the difficult noise estimation process. We can also 
mention the hierarchical decomposition | do re et al., 20 01 1 that achieves an exact power spectrum estimation 
to submaps at various resolutions, and then optimally combine them. 




(34) 
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5. Conslusions 

We have shown techniques designed to make maps from CMB data and to extract power spectra from 
them. In both cases, the brute force, maximum likelihood approach is the most correct, but generally hard 
to implement in practice. Alternative approaches, iterative or relying on Monte-Carlo simulations provide 
similar precision with smaller computer requirements. In all cases, a lot of work has to be done before : 
first by designing the instrument correctly and afterwards by cleaning the data, flagging bad samples and 
ending up with a dataset that match the minimum requirement of all the methods described in this review : 
stationarity and Gaussianity. 
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blard and P. Filliatre for reading carefully the manuscript and providing useful inputs. 
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